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We report constraints on antineutrino oscillation parameters that were obtained by using the two 
MINOS detectors to measure the 7% muon antineutrino component of the NuMI neutrino beam. In 
the Far Detector, we select 130 events in the charged-current muon antineutrino sample, compared 
to a prediction of 136. 4±11.7(stat)^g° g 2 (syst) events under the assumption |Am 2 | = 2.32xl0 -3 eV 2 , 
sin 2 (2(9) = 1.0. Assuming no oscillations occur at the Near Detector baseline, a fit to the two-flavor 
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oscillation approximation constrains |Am 2 | < 3.37 x 10 3 eV 2 at the 90% confidence level with 
sin 2 (20) = 1.0. 



PACS numbers: 14.60.Lm, 14.60.Pq, 29.27.-a, 29.30.-h 



The phenomenon of neutrino oscillations has been well 
established by experimental observations (JQ- The un- 
derlying quantum-mechanical mixing between the neu- 
trino flavor and mass eigenstates is governed by the ele- 
ments of the PMNS matrix Q , usually parameterized by 
three mixing angles and a CP-violating phase. Oscilla- 
tions are governed by the ratio of the distance traveled 
by the neutrino to its energy (L/E) and the two inde- 
pendent neutrino mass-squared differences. CPT sym- 
metry constrains the allowed differences between a par- 
ticle and its antiparticle (lpj and requires their masses 
to be identical. Differences between the measured neu- 
trino and antineutrino oscillation parameters would indi- 
cate new physics. For example, as neutrinos propagate 
through matter, nonstandard interactions fill ] could al- 
ter the disappearance probabilities of neutrinos relative 
to antineutrinos and thus the inferred oscillation param- 
Such models of new physics predict a different 
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eters 

energy dependence and so probing the standard oscilla- 
tion hypothesis to greater precision across a wide range 
of energies is valuable. 

The MINOS long-baseline experiment has made the 
most precise measurements to date of the larger (atmo- 
spheric) mass-squared splitting for both neutrinos 13 1 
and antineutrinos 0. With the NuMI facility [l5| 
configured to provide a neutrino-dominated beam, a 
measurement of v„ disappearance resulted in a mass- 
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squared splitting of Am | 
and mi xing angle sin 2 (29) > 0.90 (90% confidence limit 

[c.l.]) 0,113. 

From direct observations of disap- 
pearance, using a smaller exposure to the beam opti- 
mized for antineutrinos, MINOS measures the antineu- 
trino oscillation parameters |A77I 2 | = [3.36lQ'4o(stat) ± 

0.06(syst)] x 10" 3 eV 2 and sin 2 (20) = 0. 86±g:i2 ( stat ) ± 
0.01 (syst) (l3 |. Prior to the measurement of |A7f7 2 | by 
MINOS the strongest constraints on antineutrino oscilla- 
tion parameters came from a fit 17[ to global data dom- 
inated by Super-Kamiokande results where the sum of 
atmospheric and 77 M interaction rates was measured. 

This paper describes an analysis of the 7% V ^ compo- 
nent of the NuMI beam, optimized to produce neutrinos, 
with an exposure of 7.1 x 10 20 protons on target. The MI- 
NOS detectors are magnetized, allowing event-by-event 
separation of and 77 M charged-current (CC) events us- 
ing the curvature of the muon track. The 77^ sample 
presented here provides a new test of the oscillation hy- 
pothesis for muon antineutrinos at the atmospheric scale. 
With substantially increased statistics in the 5-15 GeV 
energy range relative to the sample obtained with the 
beam configured for antineutrinos [14| the 77^ oscillation 



probability can be probed to greater precision in this re- 
gion. 

The NuMI beam uses 120 GeV/c protons incident on 
a graphite target to produce secondary hadrons, in par- 
ticular pions and kaons of both charges. Depending on 
the sign of the applied current, two magnetic horns focus 
cither positively or negatively charged hadrons for a neu- 
trino or antineutrino beam, respectively. A 675 m long 
iron-walled decay pipe — evacuated during the first half 
of the data taking period but later filled with 0.9 atm 
helium for structural reasons — allows the hadrons and 
tertiary muons to decay in flight, producing neutrinos 
and antineutrinos. The antineutrino component of the 
neutrino beam arises from four main sources: decays of 
hadrons traveling along the axes of the horns where the 
focusing field is negligible; partially defocused hadrons 
decaying close to the horns; decays of hadrons produced 
from interactions with the helium and walls of the de- 
cay pipe; and decays of tertiary muons that arise mainly 
from decays of the focused hadrons. Muon antineutri- 
nos from neutral kaons are estimated from simulation 
to comprise 0.6% of events across the spectrum. The 
combined energy spectrum of the 77 M CC events arising 
from these sources is broadly distributed and peaks at 
approximately 8 GeV, whereas the energy spectrum re- 
sulting from the focused hadrons is narrowly-peaked at 
approximately 3 GeV. 

The two MINOS detectors [3] are located 1.04 km 
[Near Detector (ND)] and 735 km [Far Detector (FD)] 
from the target. Both detectors are segmented 
steel/scintillator tracking calorimeters. The detector 
fiducial masses are 23.7 tons and 4.2 kilotons at the ND 
and FD respectively. In CC interactions, ^(77^) + N — > 
/i~(/i + ) + X, a hadronic shower (X) and a muon track 
may be observed. The reconstructed neutrino energy is 
the sum of the reconstructed muon and hadron energies. 
Hadronic energy is measured by calorimetry. Muon en- 
ergy is measured by range for contained tracks or by cur- 
vature in a 1.4 T toroidal magnetic field for exiting tracks. 
For this data set, the fields in both detectors have been 
set so that they focus fi~ and defocus allowing the 
separation of and 77 M CC interactions. 

The inclusive 77 M CC interaction rate as a function of re- 
constructed 77 p energy is measured in each detector. The 
measured FD spectrum is compared to the projection of 
the ND data to the FD, taking into account the different 
geometric acceptances of the two detectors. In this com- 
parison, many sources of systematic uncertainty largely 
cancel due to the similarities of the two detectors. An- 
tineutrino oscillations would cause an energy-dependent 
77 M deficit at the FD compared to the projection from the 
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ND; the v ^ survival probability in the two- flavor approx- 
imation is 



Pfa -> Vp) = 1 - sin 2 (2(9) sin 2 



1.267 Am 2 L 

E 



(1) 



where L [km] is the distance from the point of antineu- 
trino production, E [GeV] the antineutrino energy, the 
antineutrino mixing angle, and Am 2 [eV 2 ] the antineu- 
trino mass-squared difference. 

Selected events must contain at least one reconstructed 
track; the longest track is identified as the muon can- 
didate. This muon candidate must originate inside the 
fiducial volume and have a positive charge determined 
from track curvature. However, the track finding algo- 
rithm can occasionally form a track out of hadronic ac- 
tivity, or misidcntify the curvature of a muon track. A 
simple charge-sign selection based on this track-fit in- 
formation yields a sample that is highly contaminated 
with both CC and neutral current (NC) events as 
shown in Fig. [TJ Monte Carlo studies show that about 
half of NC events with a reconstructed track and 7% 
of Vfj, CC events with a track are misidentified as n + 
candidates. Most of the misidentified CC events are 
high-inelasticity interactions in which the soft n~ is ob- 
scured by the hadronic shower. In addition, higher mo- 
mentum muons follow a less curved trajectory, increas- 
ing the probability of charge misidcntification. With the 
beam consisting of about 92% muon neutrinos, the ini- 
tial signal to background ratio is inherently much lower 
for muon antineutrinos than it is for neutrinos and the 
development of further selection cuts was necessary. 

To reduce the misidentified NC and CC background 
events, three selection variables are used. The first is 
a likelihood-based separation parameter based on event 
topology. The second variable is a measure of the con- 
fidence of charge-sign determination from the track fit- 
ting. The third variable provides an additional measure 
of the direction of curvature of the muon track by com- 
paring the local track direction at the vertex to that at 
the end point of the track [lj| ■ The likelihood-based sep- 
aration parameter was originally developed to distinguish 
NC background from Vu CC events in the MINOS anal- 
ysis of V/j, oscillations |l| but it is also effective in remov- 
ing the misidentified high-inelasticity CC background. 
This discriminator uses probability density functions con- 
structed from three variables: the event length, the frac- 
tion of the total event signal in the reconstructed track, 
and the average signal per plane of the reconstructed 
track. These quantities are related to the muon range, 
the event inelasticity and the average energy loss dE/dx 
of the muon track and are distributed differently for 
CC events compared to NC and misidentified CC 
events. 

The selection was optimized for statistical sensi- 
tivity to oscillation parameters equal to those measured 
for 21]. Figure [T] shows the efficiency of the full se- 
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FIG. 1: Efficiency of the selection of CC candidate events 
reconstructed with a positive charge-sign track in the Far De- 
tector. The contamination due to misidentified NC and CC 
interactions is also shown (assuming no oscillations), both be- 
fore and after all other selection criteria are applied. 



lection and the remaining contamination as a function 
of energy in the FD. Assuming no oscillations, the 
efficiency of the selection is 85% and the purity of the 
CC sample is 98%, integrated over all energies in the FD. 

The measured ND energy spectrum, shown in Fig. [2j is 
used to predict the FD spectrum, as in previous MINOS 
analyses [H[l3. 21. 



22| . This effectively mitigates sources 
of mismodcling, such as uncertainties in the neutrino flux 
or neutrino cross sections, which affect both detectors in 
similar ways. 

Hadron production in the NuMI target and beam line 
is simulated with FLUKA (H by using FLUGG as 
an interface to the GEANT4 [25J based geometry. Ad- 
ditionally, hadron production in the target is constrained 
by a fit to ND spectra [l[ , which correct the 7r and K dis- 
tributions as a function of their transverse and longitudi- 
nal momenta at production, px and p z respectively. The 
fit is performed simultaneously for several different beam 
configurations, which permits the constraint of a wide 
range of pt~Pz space for v IL parent particles. The n + /tt~ 
ratio measured by NA49 [2(| , together with the pt spec- 
tral shape from the fit, constrains the parent pt 
spectral shape, while a fit to the ND 77^ energy spectrum 
provides overall normalization and p z shape information. 
These fit parameters have been applied to the flux in ob- 
taining the simulated ND spectrum shown in Fig. [2] The 
errors obtained in the fit provide an estimate of the un- 
certainty on the hadron production from the target; the 
corresponding error on the FD event rate, extrapolated 
from ND data, is less than 1% for the beam component 
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FIG. 2: Reconstructed energy spectra of 77 M CC candidate 
events at the Near Detector. The solid line shows the Monte 
Carlo simulation, which is broken into three sources of 77 M 
parent particles. The upstream pion decay contribution orig- 
inates primarily from the target but also includes antineu- 
trinos from muons whose parents decayed in the upstream 
region. The decay pipe component corresponds to all 77 M par- 
ents (other than muons, which contribute about 3% of the ND 
spectrum) produced 45 m or more downstream of the target. 
The contribution from kaon decay is shown by the hatched 
histogram. The shaded band on the simulation shows the 
size of the systematic error on the absolute ND spectrum pre- 
diction. 



that arises directly from hadrons produced in the target. 

Figure [5] shows the contribution of different beam flux 
components to the V y, CC interaction rate in the ND 
as a function of energy. A significant fraction of ND 
events originate from parent particles produced in the 
decay pipe, predominantly from interactions of primary 
and secondary hadrons with the decay pipe walls and the 
helium (muons are not included in our decay pipe compo- 
nent definition as they are constrained by the ND CC 
events). For these events the relative acceptance of the 
ND compared to the FD is larger than for particles pro- 
duced in upstream interactions. Consequently, the con- 
tribution from decay pipe parent particles as a fraction of 
the total spectrum is larger at the ND (12%) compared 
to the FD (7%, assuming no oscillations). A systematic 
uncertainty on the size of the decay pipe component was 
assessed by scaling this component in the Monte Carlo 
simulation and comparing with the ND data. Conserva- 
tive scale factors of ±100% are applied to the decay pipe 
component, introducing an uncertainty on the total 
CC interaction rate predicted at the FD of ^50%. 

Further systematic uncertainties include a 4% relative 
normalization uncertainty between the ND and FD to ac- 
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III 


3.21 
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87 


60.2i!;£ 
76 2+ 10 9 


66.4+^ 
83.9tJJ; 6 9 


Total 


7.09 
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136.4+^;? 


150.3+^; 6 6 



TABLE I: Candidate 77 M CC events observed and expected in 
the Far Detector, broken down into two periods of approxi- 
mately equal exposure. The expected number of events in the 
oscillated case uses the parameters measured with the CC 
sample [l^ ]. 



count for uncertainties in the reconstruction efficiencies, 
exposure and fiducial masses of both detectors [2l|. A 
comparison of momentum measurement from curvature 
vs. range in stopping muon tracks constrains the uncer- 
tainty in track momentum determination from curvature 
to be 4%. The 50% uncertainty on the misidentified NC 
and CC events was estimated by scaling those com- 
ponents in the ND until the MC matched the data for 
the set of events that narrowly failed the selection on the 
likelihood-based separation parameter. The total sys- 
tematic uncertainty on the predicted number of events 
at the FD is 82% of the total statistical uncertainty, as- 
sumin g os cillation parameters equal to those measured 
for ^ 

At the FD a total of 130 selected CC candidate 
events are observed. Figure [3] shows the energy spec- 
trum of the FD data overlaid with two predicted spec- 
tra obtained from the ND data: one without oscilla- 
tions and one with oscillation parameters of |Am 2 | = 
2.32 x 10~ 3 eV 2 , sin 2 (20) = 1.0 [H. The predicted 
backgrounds arc 1.8 Vp CC events, 1.2 NC events and 
0.2 V T CC events (in the oscillated case). The integrated 
number of events observed and expected are detailed in 
Table HI The number of FD events measured in run peri- 
ods I and II is smaller than the prediction. In run period 
III, which differs due to the helium in the decay pipe, a 
larger number of events are measured compared with the 
prediction. The probability of observing a comparable or 
larger difference in event rate between the two periods, 
evaluated using mock Monte Carlo experiments, is 8.4%. 

The measured FD energy spectrum is compared to 
that predicted from the ND assuming v T oscil- 

lations, following Eq. ([T]). This comparison is made by 
minimizing a binned log-likelihood with respect to Am 2 
and sin 2 (26*). The Feldman- Cousins approach [27| is used 
to obtain confidence limits on the oscillati on p arameters 
with systematic uncertainties included [2M The con- 
fidence limits thus obtained are shown in Fig. 21 Values 
of I Am 2 1 greater than 1 eV 2 are not considered in this 
analysis, since above that point oscillations with max- 
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FIG. 3: Energy spectra of 77 M CC candidate events observed in 
the Far Detector. The predicted spectrum with no oscillations 
and with oscillation parameter values of |Am 2 | = |Am 2 | = 
2.32 x 10 _3 eV 2 , sin 2 (20) = sin 2 (26>) = 1.0 are overlaid. The 
hatched band indicates the total systematic uncertainty on 
the prediction. The estimated background includes oscilla- 
tions at the best-fit values determined by the MINOS CC 
disappearance analysis [l3|] for the CC events. 



imal mixing would cause more than 1% of the to 
disappear in the ND. Figure [4] also shows the recent 
MINOS result using the beam configured for antineutri- 
14 1, the MINOS allowed region for neutrinos fljj . 



nos 



and a fit [17( to all global data available prior to all MI- 



NOS data. The MINOS data presented in this paper 
are consistent with both the previous MINOS neutrino 
and antineutrino limits, and with the limits from a global 
fit 17j. A x 2 goodness-of-fit test using the oscillation 



parameters from [13( yields a probability of 18%. Un- 
der the assumption sin 2 (28) =1.0 these data constrain 
|Am 2 | < 3.37 x l(T 3 eV 2 (90% C.L.) in the two-flavor 
approximation. 

In summary, a high-purity sample of muon antineu- 
trino charged-current events was selected in the MINOS 
data from the 7% component of the NuMI neutrino 
beam. At the Far Detector, 130 event candidates were 
observed, which is consistent with the predicted rate in 
the case of oscillations of 136.4±11.7(stat)^g° g 2 (syst) un- 
der the assumption |Am 2 | = 2.32 x 10" 3 eV 2 , sin 2 (2#) = 



FIG. 4: Allowed regions for oscillation parameters from a 
fit to the data in Fig. [3] The region indicated by the hashing 
is excluded. Shown alongside are contours for: the MINOS 
result from the NuMI beam optimized for antineutrino 
production [l4l |; the MINOS allowed region for neutrinos [T^ |; 
and limits from a fit [l7| prior to all MINOS data. 



1.0. These data provide a new probe of the oscillation hy- 
pothesis for muon antineutrinos at the atmospheric scale. 
Significantly increased statistics in the 5-15 GeV energy 
range, compared to the sample obtained with the 
NuMI beam configured for antineutrinos, have allowed 
the oscillation probability to be measured with greater 
precision in this region and have added to constraints on 
antineutrino oscillation parameters. 
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